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ABSTRACT 

Planets with several Earth masses and a few day orbital periods have been discovered 
through radial velocity and transit surveys. Regardless of their formation mechanism, 
an important evolution issue is the efficiency of their retention in the proximity of their 
host stars. If these "super-Earths" attained their present-day orbits during or shortly 
after the T Tauri phase of their host stars, a large fraction of these planets would 
have encountered intense stellar magnetic field. Since these rocky planets have a higher 
conductivity than the atmosphere of their host stars, the magnetic flux tube connecting 
them would slip though the envelope of the host stars faster than across the planets. 
The induced electro-motive force across the planet's diameter leads to a potential drop 
which propagates along a flux tube away from the planet with an Alfven speed. The 
foot of the flux tube would sweep across the stellar surface and the potential drop 
across the field lines drives a DC current analogous to that proposed for the electro- 
dynamics of the Io-Jupiter system. The ohmic dissipation of this current produces 
potentially observable hot spots in the star envelope. It also heats the planet and leads 
to a torque which drives the planet's orbit to evolve toward both circularization and 
a state of synchronization with the spin of the star. The net effect is the damping of 
the planet's orbital eccentricity. Around slowly (or rapidly) spinning stars, this process 
also causes rocky planets with periods less than a few days to undergo orbital decay 
(or expansion/stagnation) within a few Myr. In principle, this effect can determine the 
retention efficiency of short-period hot Earths. We also estimate the ohmic dissipation 
interior to these planets and show that it can lead to severe structure evolution and 
potential loss of volatile material in them. However, these effects may be significantly 
weakened by the reconnection of the induced field. 

Subject headings: magnetohydrodynamics (MHD) - planetary systems - planets and 
satellites: dynamical evolution and stability - planets and satellites: formation - planet- 
star interactions - stars: magnetic field 
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Introduction 



A n important milestone in planetary astronomy is the discovery of a Jupiter-mass planet, 51 

Peg b (jMavor and Queloz 19951 ). Its extraordinary 4 d ay orbital period rekindled a theoretical ex- 

pectation that protoplanets may undergo orbital decay (jGoldreich &: Tremaine 198Q ; iLin &; Papaloizou 19861 ) 
as a consequence of their tidal interaction with their natal disks. Today, more than 500 planets 
have been discovered around nearby stars. Among them, there is a pile up of ~100 planets with 
periods (P) less than a week and mass (M p ) two orders of magnitude larger than that of the Earth 
(-M©). Transit observations of some of these planets indicate that they have a radius and density 
comparable to that of Jupiter and Saturn and are commonly referred to as hot Jupiters. 

These hot Jupiters, including 51 Peg b, w ere formed in pre s umably at some preferred locations 
(beyond the snow line) of their natal disks (jlda fc Lin 20081 ; llda fc Lin 20101 ). After acquiring 
sufficient masses to open gaps in their natal disk, they undergo type II migration along with the 
viscous diffusion of their surrounding gas. In order to account for their survival, we proposed that 
the migration of 51 Peg b and other short-period gas giant planets may have s talled when the y 
entered into the magnetospheric cavity of their host star during their infancy (ILin et al. 19 96). 
Since protostellar disks are expected to be truncated within the magnetosphere of their central 
stars (see below), once any protoplanet enters into this region, its migration would slow to a halt 
as the intensity of its tidal interaction with its nascent disk weakens. 



The existence of magnetosphere around T Tauri stars was pr oposed (IKonigl 199lh to account 
for the observed period distribution which peaks around 8 days ( Bouvier et al. 19931 ). If efficient 
angular momentum flow between protostellar disks and the magnetosphere of their central stars 
can enforce corq tation at their interface where the magnetic and viscous torques are balanced 
(jShu et al. 19941 ). we would infer kilogauss fields. Today, Zeeman splitting of emission lines have 
been directly measured for many T Tauri stars an d these observations confirm the presence of 
kilogauss fields on their surface (|Johns-Krull 20071 ). If this field strength represents that of a 
dipole stellar field, the magnetospheric radius for T Tauri stars with accretion rates in the range 
of lO~ 8 folO -7 M0 yr -1 would extend to regions beyond the orbits of their close-in planets. The 
radius of this magnetospheric cavity expands during the depletion of the disk gas and the decline 
of the accretion flux through the disk. 

Although the magnetospheric-cavity scenario provides a useful qualitative model for the abun- 
dant population of short-period planets, a detailed reconstruction of the observed period distribu- 
tion requires a quantitative determination of short-period planets' retention efficiency. After they 
have entered the magnetospheric cavity of their host stars or after they have been engulfed by the 
magnetospheric cavity, protoplanets continue to interact with the stellar magnetic field. It is not 
clear whether this process can induce planets to undergo further orbital interaction. 

In order to und erstand the nature of this physical mechanism, we carry out a series of investi- 
gations. In Paper I (jLaine et al. 20081 ). we provided a general description of the relevant physical 
effects. We first consider the retention of short-period gas giant planets inside the stellar magne- 
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tosphere. We decompose the stellar magnetic field into a steady and a periodically modulating 
component and the planet into a day and night side. In this previous investigation, we considered 
the periodically modulating component of the field which can be due to either the planet's eccentric 
orbit or the star's non-synchronous (with respect to the planet's orbital angular frequency) and 
non-aligned (with respect to the star's magnetic poles) spin. On the night side of the planet where 
the magnetic diffusivity is relatively high, this time-dependent field can permeate into the planet's 
envelope and induce an AC current. The ohmic dissipation of this current not only heats the 
planet but also provides a torque which drives the planet's orbit toward a state of circularization 
and synchronization with the star's spin. We also found that close to the host star where the stellar 
magnetosphere is intense, ohmic dissipation can cause a planet's interior to heat up such that it 
expands and overflows its Roche lobe. The gas flow from the planet to its host star via th e inner 
Lagrangian point also transfers angular momentum to the orbit of the planet (|Gu et al. 20031 ). This 
process can halt the migration of the planet despite its loss of angular momentum as a consequence 
of its tidal interaction with the disk and the host star and the direct torque applied by the stellar 
magnetosphere on it. The results of this analysis will be presented elsewhere. 

Another class of planets have been discovered with M p ~ a few and P in the range of 
a few days to 2 months. In contrast to their Jupiter-mass siblings, these pla nets probably have 
rocky or icy interna l structures and are commonly referred to as super-Earths (jMavor et al. 20081 ; 
Howard et al. 20101 ) . Recently Kepler mission led to the discovery of over 1200 planetary candi- 
dates. Since they all have radii more than twice that of the Earth, they may also be super-Earths, 
albeit their masses are yet to be determined. These super-Earth s too are probab l y formed at loca - 
tions ranging from a fraction to several AUs from their host stars (jlda Lin 20081 : llda Lin 201Q ) . 
In contrast to the gas giants, super-Earths may not have adequate mass to open a gap and undergo 
type II migration. Nevertheless, they do tidally interact with their natal disk and undergo type 
I migration, due to an imbalance between the torque they exer t on the disk at the Lindblad and 
corotation resonances both interior and exterior to their orbits (jTanaka et al. 20021 ). 



In most regions of the disk, type I migration is directed inward. However, due to the corotation 
torque, there are m i gration barriers where the orbital decay of isolated super-Earths may be halted 
(jMasset et al 20061 : iPaardekooper et al. 20101 ). These barriers occur at gas surface density (S g ) 
maximum near the snow line (a; ce where water vapor condensate), at the inner edge of a dead 
zone (adead where the ionization fraction near the mid plane is too small to maintain coupling 
between turbulent magnetic field and d isk gas), and just outside the magnetospheric cavity, a mag 
(|Kretke fc Lin 20071 : iKretke et al. 20091 ). 



Since E g and the mid-plane temperature of the disk decline with time, the location of these 
barriers also evolves with them. During the advanced stages of disk evolution, S ff may be sufficiently 
small that the dead zone essentially vanishes and the stalled embryos near ai ce and adead resume 
their migration. Assuming the field strength decreases slowly or remains unchanged, the location of 
o mag also expands beyond the orbital semimajor axis of the super-Earth s, ose, during the tra nsition 
from classical to weak-line T Tauri phases on a timescale of ~ 10 Myr (jKretke fc Lin 20101 ). After 
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the disk dep letion, the stellar magn etic field and spin rate decrease (jSkumanich 19721 ) on a timescale 



~ 100 Myr (jSoderblom et al. 19931 ). During these evolutionary stages, the location (a cor ote) where 
the frequency of the Keplerian motion matches that of the stellar spin also evolve relative to ose- 
Thus, we anticipate that some super-Earths are located interior to the corotation radius while 
others are located beyond it. 

In this paper, we consider the interaction between close-in super-Earths with the steady com- 
ponent of their host stars' magnetic field. Beside differences between their masses and radii, rocky 
planets have much larger electric conductivity on their surface than the envelope and atmosphere 
on the night side of close-in young gas giant planets. The dayside of gas giant planets is exposed to 
the stellar ionizing photons and may have much higher ionization fraction and electrical conductiv- 
ities than their night side. We will consider this more complex aspect of the hot Jupiter problem 
elsewhere. 

The super-Earth problem we are considering here is analogous to the Echo satellite (in the form 
of a conductor) moving relative to background (Earth) field line which was analyzed by Drell et al. 
(1965). In general, an electric field is induced across the conductor (in directions orthogonal to the 
field and motion). However, the electric field must vanish (in its frame) inside a perfect conductor 
(with an infinite conductivity). On its surface, the electric field generate a current which leads to 
an induced magnetic field. The induced field cancels the unperturbed field inside the conductor 
so that there is no relative motion between the perfect conductor and the net (unperturbed plus 
induced) field inside it. Outside the moving conductor, the net field appears to wrap around it. 
At large distances from the moving conductor, the induced field propagates away from it with the 
Alfven speed. 

In a slightly different context, Goldreich and Lynden-Bell (1969) analyzed the electrodynamic 
interaction between Io and Jupiter. They treated Io as a conductor moving in Jupiter's magneto- 
sphere. Because Io's conductivity is larger than that of Jupiter, it drags a flux tube of field lines. 
They showed that an electromagnetic field is induced across Io's surface. Provided the conductivity 
is high along and low across the field lines, the electric potential drop across them would propagate 
and be maintained along the field lines connecting Io and Jupiter with an Alfven speed. At the 
foot of the flux tube where it enters into Jupiter's atmosphere and envelope, conductivity across 
the field lines increases with the density of the surrounding gas such that the potential drop drives 
a DC current across the potential drop. Provided that the induced field can propagate back to Io 
before the unperturbed field slips through it, the current forms a close circuit. In this scenario, Io 
acts unipolar inductor. 

In this paper, we apply the Goldreich and Lynden-Bell's model to the study of a super-Earth 
moving in its host star's magnetosphere. For computational simplicity, we neglect the planet's 
intrinsic magnetic field, as we have done in the previous paper. In the problem we are considering, 
a steady component of the stellar magnetic field is present regardless of the differential motion 
between the star's spin and the planet's orbit. In an asynchronous system, a flux tube of magnetic 



- 5 - 



field between the planet and the star cannot be infinitely anchored on both entities. If the planet's 
conductivity is smaller than that on the surface of the host star, the flux tube would tend to be 
anchored on the surface of the star along with all other unperturbed field lines and it would slip 
through the planet. 

In Section 2, we introduce qualitative discussion and schematic illustration to show in the limit 
that the electric conductivity in the super-Earth planet is higher (but not by an infinite amount) 
than that of its host star's envelope, the relative motion between the planet and the stellar magnetic 
field leads to an induced emf and a potential drop across the planet. Outside the planet, a flux 
tube of (unperturbed plus induced) magnetic fields would appear to be approximately anchored on 
the planet. In the tenuous regions between the planet and its host star, the conductivity along the 
field lines is much higher than that across them and the electric current flows freely to maintain 
constant electric potential along them. In the absence of field reconnection, Alfven waves propagate 
to infinity on open lines and the electric current flows to the surface of the host stars along closed 
field lines. Due to the finite resistance of the surrounding (stellar) gas, the foot of the flux tube 
on the surface of the host star slips through the stellar atmosphere and the electrical potential 
drop across the foot of the flux tube drives an electric current with an associated rate of ohmic 
dissipation. 

Based on the above qualitative scenario, we construct a quantitative model in this paper. In 
Section 3, we introduce the values of the different parameters we use in the numerical applications 
and derive the analytical expressions for the induced electric field, intensity, ohmic dissipation, and 
torque. In our numerical applications, we adopt a set of fiducial physical parameters for a rocky 
super-Earth planet with two Earth radii, 2i? ffi ps 1.4 x 10 7 m on a 3 day circular orbit (a ps 6 x 10 9 m 
around a T Tauri star with a mass equals to that of the Sun (1M ps 2 x 10 30 kg) and a radius twice 
that of the Sun R* = 2R Q ps 1.4 x 10 9 m. We also assume that it has a solar luminosity (ILq), a 
surface temperature of T* = 4 x 10 3 K, and a dipole field with a strength 0.2 T (i.e. 2 x 10 3 G) on 
the stellar surface (which corresponds to a magnetic dipole strength of 5.4 x 10 33 Am 2 ). 

With these parameters, we first analyze two limiting cases of rapid and slow stellar spin. We 
discuss the condition of validity of the model in Section 4 and derive in Section 5 the expressions and 
values of the resistances and Alfven speed. We will consider more general sets of model parameters 
elsewhere. In Section 6, we calculate the values of the induced intensity, ohmic dissipation, and 
torque, and discuss the relevance of these values. In the context of planetary migration in the 
presence of their natal disk, we show that if the planet orbits around the host star outside its 
corotation radius (i.e. the Keplerian frequency of the planet's orbit is smaller than the stellar spin 
frequency), the net torque associated with this induced current would provide an adequate rate 
of angular momentum transfer to balance against the rate of tidally induced angular momentum 
loss by the rocky planet to the disk. In the limit that the planet is inside its host star's corotation 
radius, the planet's orbit would continue to decay. Finally, we summarize our results and discuss 
their implications in Section 7. 
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2. Qualitative Illustration and Description of the Phenomena and Brief Outline of 

the Calculation 

We consider a rocky planet with a mass of several M e moving in the dipole magnetic field of 
the star it is orbiting. The relative motion of such a conductor in an external stellar magnetic field 
generates an induced emf across the planet. There are two complementary effects (see paper 1) 
described by the complete magnetohydrodynamic (MHD) equation: the diffusion of the magnetic 
field in the planet and the magnetic induction (with its associated drag). 

The MHD equation describing the electrodynamics of the planet in the stellar field can be 
written as 

f)T3 

- = VA(uAB)-VA(ijVAB) (1) 
where rj is the magnetic diffusivity (which is equal to (fi^a)^ 1 where a is the electric conductivity). 

In Paper 1, we focused on the diffusion of the stellar magnetic field inside a hot Jupiter. In 
order to do so, we considered the periodic component of the stellar magnetic field felt by a planet 
when the axis of the stellar magnetic moment is not aligned with the planetary orbital axis, or when 
the planet is on an eccentric orbit. The drag of the magnetic field by the planet due to induction 
can be neglected if the planet's orbit co-rotates with the star's spin or if the electrical conductivity 
of the planet is low compared with that of the outer layers of the star (which is the case at least in 
the night side of a hot- Jupiter). 

In this limit, the diffusion of the stellar magnetic field inside the planet is modulated by the 
electric conductivity (inversely proportional to the magnetic diffusivity) profile in the planet. A 
somewhat higher electric conductivity in the planet would tend to decrease the penetration depth 
of the stellar magnetic field and the volume where the electric current induced by the field can be 
dissipated, but it would also increase the volumic ohmic (power) dissipation rate. Likewise, a lower 
conductivity in the planet would enable the stellar magnetic field to penetrate deeper into it, albeit 
the induced current also encounters a lower volumic ohmic dissipation rate. Consequently, we found 
that, the total ohmic dissipation rate over the entire planet does not change significantly over a 
reasonable range of electric conductivity for a hot-Jupiter (neglecting the effect of photo-ionization 
in its atmosphere). 

In the present paper, we describe the induction (and associated "drag" of the field lines) which 
was neglected in paper I. For simplicity, we neglect the damping of the magnetic field in the planet 
associated with the diffusion term and focus on the case where the planet is able to significantly 
drag the stellar magnetic field lines which are enclosed in the flux tube that passes through the 
planet. 

Throughout this paper, this "flux tube which passes through the planet" is simply referred to 
as "the flux tube" (we are interested in the part that extends between the interior of the planet and 
above the surface of the star). The "foot of the flux tube" refers to that part of the flux tube which 
extends below the surface of the star for a distance D pn (the subscript " pn" refers to penetration) 
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to be estimated below. 

Between the planet and the surface of the star, the volumic current flows along the flux tube 
(parallel to the magnetic field lines-the electrons in fact gyrate around the magnetic field lines) , but 
the volumic current crosses the flux tube in the planet and at the foot of the flux tube (perpendicular 
to the magnetic field lines in the stellar atmosphere). Figures 1 and 2 present the general overview 
of the system. 

In this model, we make a distinction between the regular (unperturbed) stellar dipole magnetic 
field (magnetosphere) which co-rotates with the stellar spin and the field lines which define the flux 
tube (composed of both the stellar field in the planet and the induced flux tube) which appear to 
be dragged along with the planet and thus move relative to the rest of the stellar magnetosphere. 
This drag is significant whe n the electric conductivity of the planet is large compared to that of 



the outer layers of the star (jPiddington Sz Drake 19681 ) 



A large conductivity in the planet's interior would lead to a large surface current which induces 
a field and cancels the external (unperturbed stellar dipole) field. With small magnetic diffusivity, 
the planet's interior would appear to be shielded from any time dependent external magnetic field. 
We, therefore, consider only the time-independent component of the stellar magnetic field. In this 
model, the time- independent component of the stellar magnetic field permeates the planet. The 
motion of the planet relative to the stellar magnetosphere induces an electric field £, an induced 
volumic electric current S, and an induced difference of potential U across the planet's diameter. 
In Figure O we provide a schematic illustration on the field lines and current inside the planet. 

For the flux tube between the planet and its host star, we adopt the assumption of high electric 
conductivity along the magnetic field lines and low conductivity across them which was introduced 
by Goldreich & LyndenBell (1969) for the system Io-Jupiter. In this case, the difference of electric 
potential (induced by the planet's relative motion with respect to the magnetic field of its host star) 
drives an electrical current out of the planet, along the flux tube, across its foot on the atmosphere 
of the star, and back to the planet along its other half (Figure 1). 

In the limit of negligible electrical conductivity in the direction normal to the flux tube, the 
electric current can only cross the field lines in the planet or in the atmosphere of the host star. 
The assumption of high electric conductivity along the magnetic field lines also implies 1) that the 
difference of potential U across the planet's diameter is transmitted along the flux tube without 
significant drop in potential and 2) that the plasma enclosed in the flux tube is dragged along by 
the motion of the magnetic field lines. An electric circuit is therefore created, where the flux tube 
acts as electric wires, the planet as a unipolar inductor with internal resistance 1Z P , and the foot 
of the flux tube as the largest resistance. 

We see that there are in fact two circuits (see Figure 1). The first one is composed of the foot of 
the flux tube on the northern hemisphere of the star, its corresponding flux tube, and the northern 
hemisphere of the planet. The second is equivalent and symmetric to the first one (the plane of 
symmetry being the plane of the planetary disk). Except when explicitly stated, the calculations 



- 8 - 




Fig. 1. — Schematic illustration of the planet-magnetosphere interaction system. The circuit di- 
agram idealizes the basic physics which is described in the text. In this context, a planet with a 
non-negligible motion (into the plane of the diagram) relative to the stellar magnetosphere induces 
an emf. At the location of the planet, the direction of the unperturbed stellar dipole field is point- 
ing downward. The potential difference across the flux tube generate a current with a flux which 
is primarily determined by the electrical resistivity in the atmosphere of the host star. Arrows 
indicate the flow direction of the current. 



- 9 - 




Fig. 2. — Main parameters in the calculation of y\ and yi which defines the geometry of the foot 
of the flux tube 
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Fig. 3. — External and induced field around a planet moving relative to the magnetosphere of its 
host star. The planet is moving out of the plane and toward the lower left side of the illustration. 
Arrows in this idealized cartoon illustrate the direction of current flow across the potential drop 
due to the induced emf. The current also lead to an induced field which distorts the field near the 
planet. Information on the induced emf propagates along the field lines in the direction of the host 
star. 
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(current, resistances, ohmic dissipation, torques, etc.) describe only one circuit (the northern one). 

In an electric circuit composed of a generator (with an emf j Sdl and resistance 1Z 9 over a 
length scale I) and other resistances along the circuit 1Z (here primarily the resistance of the foot 
of the flux tube), the intensity of the current / is determined by J Sdl — !Z g I j Sdl = U = 1ZI. 
With the parameters we have adopted here, we show in Section 5 that the resistance along the flux 
tube T^tube an d the resistance across the planet 1Z P are small compared to that across the foot of 
the flux tube on the star 1Z*. In this limit, 1) the potential drop across the planet with a radius R p 
is U ~ 2£R P , 2) the magnitude U is approximately constant along each field lines in the flux tube 
between the planet and its host star because the resistance of the tube 7£ tu be an d the induction 
are negligible, and 3) the total current is determined by the largest resistance along the circuit, i.e. 
that at the foot of the flux tube in the stellar atmosphere. 



3. Description of the Model: Values of the Parameters, Geometry, Analytical 
Expressions, and Equation of Evolution of the Planet's Orbit 

3.1. Values of the Parameters for a Fiducial Model 

Except when explicitly stated otherwise, we consider a system composed of a super-Earth 
orbiting closely a young T-Tauri star with a time-independent magnetic dipole. We assume that 
the magnetosphere co-rotates with the star, and adopt the following numerical values (SI units) for 
the parameters intervening in the model: 

For the star we adopt the following model parameters: 
Temperature of the isothermal outer layer: T* = 4000 K, 
Radius: R+ « 2R & 1.4 x 10 9 m, 
Mass: M* ~ M ~ 2 x 10 30 kg, 

Opacity at the photosphere: k ~ 3 m 2 kg -1 (which is equivalent to taking a surface pressure of 
about 15 Pa), 

Magnetic dipole strength: m = 5.4 x 10 33 Am 2 , w hich corresponds to a magnetic field of 0.2 T 



(Tesla) = 2 x 10 3 Gs (Gauss) at the stellar surface (jYang et al. 20081 ). and 
Spin period: 8 days (slow-rotator) or 0.8 days (fast-rotator). 
We will consider more general stellar models elsewhere. 

For a super-Earth, we consider the following case: 
Radius: R p ~ 2i? e ~ 1.4 x 10 7 m, and 

Semi-major axis: a ~ 0.04 AU ~6x 10 9 m (which corresponds to a period of 3 days). 
In the electrodynamics of super-Earths, the magnitude of M p does not enter explicitly (it does 
implicitly through the radius) the calculation of the torque, albeit their orbital evolution timescale 
does depend on it (for example, see Equation (j!6p . The linear speed in a frame co-rotating with 
the star (v p / s = (uj p — w*)a) of such a planet orbiting a star rotating slowly is 9 x 10 4 ms _1 and 
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4 x 10 5 ms 1 around a star rotating fast (these are the absolute values). 

3.2. Length and Width of the Foot of the Flux Tube in a Spherical Approximation 

As defined above, the "flux tube" refers to the flux tube composed of the field lines of the 
magnetosphere that pass through and are dragged along by the planet. This flux tube connects 
the planet and the surface of the star, and its foot penetrates into the star to a depth which will 
be determined later. 

The stellar magnetic field has the geometrical structure of a magnetic dipole field. Thus, the 
foot of the flux tube at the stellar envelope (the circled area in figure [2]) is, at the first order in 
y/R*/a, an ellipse which axes have respective lengths 

(R \ 3/2 

y 2 = 2i? P MJ , (3) 

sinOp = \ — (4) 

V a 

cosOp = \1 = s, (5) 

V a 

where Bp is the angle between the stellar spin axis and the location of the foot of the flux tube. 
When the current J crosses the foot of the flux tube in the stellar envelope, it covers a length y\. 
In the rest of the paper, we take cosOf is roughly equal to 1. We represent it with the symbol s in 
analytical equations and take it to be equal to 1 in numerical applications. 

In Figure HI we zoom in on the foot of the flux tube at the stellar atmosphere. In order to 
derive y\ and y 2 (at the first order in R*/a), we first solve B A dl = and obtain sina = a +R p 

and sinj3 = \J a ^* R with a and j3 defined in Figure El We then write y\ = R*(a — f$) and 
2/2 = |^27r R^sinO 'p. 

For a super-Earth (using the model parameters listed above), we find y± = 1.6 x 10 6 m and 
?/2 = 3.2 x 10 6 m. For a hot-Jupiter, we would typically need to multiply these values by a factor 10. 
The height of the foot of the flux tube is derived below in Section T5. 3. 31 The numerical applications 
are for s = 1. For semi- major axes comparable to the stellar radius, the multiplicative factor 
(1 — R*/a)~ l l 2 in y\ would significantly affect the length of the foot of the flux tube. 



Fig. 4. — Penetration of the flux tube in the stellar atmosphere. Potential difference across the 
planet would be maintained at its foot print on the surface of its host star if there is sufficient time 
for the Alfven waves to transit this information. Arrows indicate the flow direction of the electrical 
current. The values of y\ and 2/2 are given by Equations The top circle represents the stellar 
surface and the penetration depth D pn = — r pn is evaluated in Section 15.31 
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3.3. Induced Difference of Potential 

The planet is a conductor moving in the stellar magnetosphere with relative linear speed v p i s . 
Modeling the stellar magnetic field as the one created by a magnetic dipole (of magnetic moment 
m), the magnitude of the induced electric field in the planet £ p is 

£ P = v p / s B*(a) = (u} p - L0 *) a ^^ ( 6 ) 

with jB*(a) represents the stellar magnetic field at the location of the planet. Numerical applications 
for a super Earth give £ p = 240 V m _1 (slow-rotator), 1000 Vm _1 (fast rotator) (a slow or fast 
rotator depends on the spin frequency of the star, as described in section 3.1). 



The magnitude of the difference of potential U (or emf ) generated across the planet is thus 

U = 2Rp£ p = 2Rp(u> p — w*)a— — (7) 
' ' Ana 6 

For the super-Earth models under consideration, U = 6.7 x 10 9 y (slow rotator) and 2.8 x 10 10 y 
(fast-rotator). This difference of potential is transmitted across the flux tube (with the assumption 
of infinite conductivity along the flux tube that passes through the planet) and generates a uniform 
electric field £* in the stellar envelope (see Figure 1) at the foot of the flux tube 

U nom ( a \ 3/2 

The value of £ * does not depend on the radius of the planet, and for our values of the parameters 
for a young T Tauri star, £* = 4.2 x 10 3 Vm (slow rotator), 1.8 x 10 4 Vm _1 (fast rotator). 



3.4. Analytical Expressions: Intensity in the Circuit, Ohmic Dissipation, and 

Torque in the Planet and Star 

The induced current I is given by 

I = Jdydz = £*y 2 / a*(z)dz = U— / a*(z)dz = 4i?p(w p - ui^a ^^ s / a*(z)dz. (9) 

In the previous equation, J is the volumic electric current in the stellar atmosphere at the foot of 
the flux tube (induced by U), y varies from to yi (the width of the foot of the flux tube), and z 
varies from r pn (the radius to which the flux tube can penetrate into the stellar atmosphere) to 
and H = 2s (see Figure In this circuit, the total resistance is the sum of that across the planet, 
the foot print of the flux tube on the stellar surface, and along the flux tube. Here we consider 
only the largest contribution and neglect that across the planet. In Section [5.3^ we determine the 
magnitude of r pn and evaluate 

a p (z)dz (10) 
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such that I = 2UT,s. In the above equation, <J p {z) is the local Pedersen conductivity (see Equation 
|4"T|) . The total resistance of the stellar atmosphere at the foot of the flux tube is 

U yi 1 1 
= — = — — = . (11) 

The total ohmic power dissipation in the stellar atmosphere at the foot of both flux tubes (one 
for each hemisphere, thus the multiplicative factor 2) V* and in the planet V p are 

= 21ZJ 2 = 2UI = 4C/ 2 £s = 16i^( Wp _^)V (|^) 2 £s (12) 
V p = 2U P I 2 = 21Z P {2U^sf = vj^- (13) 

where 1Z* is the resistance of the foot of the flux tube (on one hemisphere), and 1Z P is the resistance 
of one hemisphere of the planet. Since the resistance in the planet is considerably smaller than 
that in the star, most of the power is dissipated in near the foot of the flux tube on the surface of 
the star. Note that the magnitude of both V* and V p is determined by I and 1Z* . 

The total torque (for both circuits, one circuit for each hemisphere) due to the Lorentz force 
(the axis is the stellar spin axis) on the star (equal in absolute value to that on the planet) is 

% = 2 jf J r A (J A B)dxdydz = 2(R* S ine F )(yiI)( 2 ^ ™™f F ) = AR p a^Ie z , (14) 

% = im 2 p a 2 (u p -u*) (£^) 2 £se z (15) 

where s = cosQf as defined in ([2} and I is the integral of the volumic current across a cross section 
(we take an averaged view of the volumic current rather than determining its complex geometry 
inside the planet). We have calculated here the total ohmic dissipation and torque (i.e. for both 
hemispheres). 



3.5. Equation of Evolution of the Stellar Spin and Planet's Orbital Angular 

Velocity and Semimajor Axis 

The torque on the planet T p is equal and opposite that on the star %. Consequently, the 
semimajor axis of a super-Earth on a circular orbit evolves at a rate 

2a ^ 

a = —T p (16) 
tip 

where the total angular momentum of the planet's orbit is H p = M p a 2 ujp. Since the total angular 
momentum of the system is conserved, the changing rate of the stellar spin is 

= zkm (17) 
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where c* ~ 2/5 is the inertial constant of the star. According to the above expression, the planet 
would undergo orbital decay and its host star would spin up if it is inside corotation (or equivalently 
if ojp > a;*). Similarly, the planet's orbit would expand and its host star would spin down if it is 
outside corotation. 

The planet's orbital frequency uj p is related to the semi-major axis oj p = ^/sf * ■ Using the 
expressions calculated for the torques and Equations (|16p and (|17p . we find 



■ _ M P (GM^ oj p 10 uj p 



OJ 



where the numerical application was for a planet with 10 Earth Masses (M p = 6 x 10 25 kg) and 
s represent second. Therefore, we can estimate the variation of during the evolution of the 
planet's migration 

. . M„(GM») 2 / 3 . A / 1 \ , 

i^^-m-^W 1 <i9> 

which is negligibly small. 

We can thus consider that the star's angular velocity is roughly unaffected by this transfer of 
angular momentum. Using Equation (|16[) . a 3 = ffi , ^ = -jf-^f-, and the expression for the torque 



on the planet (equal to —% with % given in Equation (fT5|) ) we find 



^m 2 ^=a (20) 



oj*(oj p - w*) M p V 4vr / (GM 



expression for a is defined as in[2l The previous equation becomes 



where oj* is constant and a = 48 (^j?) (gfjfp (°f um t s 3 ), and s in the numerator of the 



1 ojp 1 oj p 1 ojp 1 cJp 1 



OJ* UJ. 



1 ojI ojj uj'I wj w £ w p w| ( Wp _ w 



• -— ^— T «• (21) 



p ^* ^p 



After integration, we find 



—^ + —^ + m+ln(\l-m |) = -J. + -J. + n i + i„(| l-fl, \) + auji{t-t ) (22) 

where we define D,(t) = uj*/oj p (t) and fij = 0(t = to). 

Near co-rotation (i.e. f2 ~ 1), the In function is dominant and we find 

ln(\l- tt(t) |) = In (| 1 - a |) + aujf(t - to). (23) 
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If oj p is greater than , then 



UJp(t) 



fi(t) 



I" 1 



,*. s ) exp (aui(t - t )) 



(24) 



Similarly, if u) p is less than w*, then 



n(t) 



1 + 



Up (to) 



l) exp(a^(t-t )) 



(25) 



and the planet-star system approaches co-rotation with a timescale 



1 



M P (GM*) 2 f^om 



-2 



In the general case, 0J p {i) follows Equation (|22p . or written differently: 



(26) 



f(t) -ln\l 



UJ p (t) 



<(t)-^p(t) = ^ 



where /(t) = A + 5(t - t ) with A = + H?/2 + n< + In (| 1 - n< |) and 5 



aw; 



(27) 



In order to get an equation of evolution of the semi-major axis, one can replace in Equation 
Z3l near co-rotation) or ([22l in the general case) O p (t) by (a(t)/a c ) 3 ^ 2 and by (a(to)/a c ) 3 ^ 2 with 
a r the co-rotation radius. Near co-rotation, we find 



a(t) = a c 



a(to' 



3/2' 



exp 



t - to 



2/3 



(28) 



for the case where a is smaller than a r 



4. Condition for the validity of the Model: t^ < i maa; 

In order to apply the model described above, one needs to verify that the time tA required for 
the Alfven waves to travel along the flux tube (to a depth D pn inside the star to be determined), 
and back to the planet is smaller than the time t max it takes the flux tube to slip ahead of the 
planet by more than its diameter. This condition ensures that a perturbation along a field line of 
the flux tube has the time to travel back and forth while the field line is still part of the flux tube 
that passes through the planet. Figures [5] and [6] illustrate this condition. Figure [7] shows the field 
lines near the planet in the case where the condition is not met. 

In order to calculate tA, we need to estimate the Alfven speeds along the flux tube between 
the planet and the star and in the stellar ionosphere (at the foot of the flux tube). Similarly, the 
calculation of t max requires the value of the conductivities (or resistance) of the different components 
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Scenario 1 : the planet is too far away. t A , tube>t 




Fig. 5. — Propagation of the Alfven wave between a relatively distant planet and its host star. 
The planet's motion relative to the stellar field induces a potential drop across the flux tube in the 
proximity of the planet. This information propagates along the flux tube toward the host star with 
an Alfven speed. Due to finite diffusion and the relative motion between the planet and the stellar 
magnetosphere, the net field lines also slip through the planet. In this illustration the timescale 
required for the Alfven wave to reach the host star is long compared with that required for the 
slippage of the field. The circuit is not established in this case. 



Scenario 2: the planet is sufficiently close 

Planet 




Fig. 6. — Necessary condition for a complete unipolar inductor circuit. Similar to the illustration 
in Figure [5l the potential drop across the planet propagates along the flux tube toward the planet 
with an Alfven speed. In this illustration, the planet is sufficiently close to its host star that the 
potential drop can be established on the surface of the host star before the fields slip through the 
planet. This potential drop induces a current which is determined by the resistivity on the stellar 
surface. In this case, it is possible to complete the circuit induced by the motion of the planet. 
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Fig. 7. — Magnetic field lines in the case where the circuit is not closed. The planetary motion 
relative to the stellar magnetic field induces charge separation. Without any connection (at infinity) 
between the separated charges, there is no current in the frame of the moving planet. In a infinitely 
conducting planet, the separated charges are concentrated near its surface. The magnitude of the 
induced field £ p is determined by Equation ©. In the stationary frame (centered on the host star), 
the separated surface charges carried by the planets generate two opposite currents as well as a 
finite V A£ p in the moving planet. Interior to the infinitely conducting moving planet, the induced 
field exactly cancels the unperturbed field as if there is no diffusion of the stellar magnetic field into 
the planet. In the external region close to the moving planet, the induced field strongly perturbs the 
stellar magnetic field. The net field distortion is symmetric around an infinitely conducting moving 
planet (analogous to an invisic flow around a spherical object) such that there is no net torque 
acting between the planet and the stellar field. This symmetry would be broken and the drag on 
the planet would be finite if its conductivity is sufficiently low to permit significant slippage of the 
stellar magnetic field or if a complete circuit connecting two sides of the planet can be established 
on the surface of its host star. 
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of the circuit. Indeed, the ratio of 1Z P (resistance of the planet) and 1Z* (resistance of the foot of 
the flux tube in the stell ar atmosphere) determines the amount of relative slippage between the 



flux tube and the planet (jDermott 19701 ) . 



In the limit where 7Z p is comparable to or larger than 7Z* (as in the night side of synchronously 
spinning hot Jupiter, see Paper I), the flux tube would tend to slip through the planet. In this 
case, the flux tube would slip ahead of the planet by a distance ~ 2R P in a relatively short time 
tA, and it might not be possible to maintain a closed circuit. 

In the most unfavorable case {7Z P /TZ* = oo, i.e. the flux tube passes through the planet 
completely undisturbed), tA = 2R p /v p / s (where v p / s represents the speed of the planet in the 
frame rotating with the stellar magnetosphere) . In the opposite extreme limit, 1Z P /1Z* = and 
tA = oo such that the flux tube is completely anchored in the planet. Differential motion steadily 
stretch the field lines until they reconnect. A more realistic situation falls somewhere between these 
two extreme limits, and the smaller 1Z P /1Z*, the easier it is to satisfy the condition of validity. 



4.1. Qualitative Estimate of the Relative Slippage 

A first qualitative criterion is given by the following argument. We want to determine whether 
the magnetic flux tube slips on the planet or on the star. In the absence of any companions, the 
magnetic field in the magnetosphere of a star rotates with the star. A close-in planet would tend 
to drag the stellar magnetic field lines that pass through it along with its motion. 

Let us define u p , w*, and ujb to be the angular velocity of the planet, star, and magnetic field 
in an absolute frame. We then consider the relative motion between the planet and the field lines, 
and between the field lines and the star. We write £l p = uj p — ujb and il* = ujb — and our 
goal is to estimate the magnitude of fip/fi*. In these notations, the planet's speed relative to the 
magnetic field is then v p = £l p a (a is the semi-major axis), and the speed of the field lines (that 
pass through the planet) relative to the star is v* = Q*R* (R* being the stellar radius). 

Considering the DC component of the field, we can write the complete MHD induction equation 
(see Equation ([1])) for the star and the planet. 



dt 
dB* 



V A (v p A Bp) + V A (??pVA B % 



V A (u* A B*) + V A (??*VA B*). (29) 



m 

In a steady state, the first and second equations imply v p w r] p /R p and v* w r]*/R*, respectively. 
Therefore, we obtain 



(30) 



In the context of Io- Jupiter interaction R^ /R p a ~ 7. Since the electrical conductivity on Io is 
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estimated to be much larger than that on Jupit er, the flux tube which passes th rough Io moves with 
Io and drags its foot on the surface of Jupiter (jGoldreich &: Lvnden-Bell 19691 ) . For a Jupiter-mass 
planet orbiting a young T-Tauri star (with radius twice that of the sun) at 0.04 AU (3 day-period), 
Rl/Rpd ~ 4, and the anchorage of the flux tube which passes through the hot Jupiters is thus 
determined by the ratio of the diffusivity through the planet to that through its host star. For a 
super-Earth with radius twice that of the earth at 0.04 AU, R\j 1 R p a ~ 20. 



4.2. Analytical Expression of the Time Constraint, t max 

The magnitude of t max is the time it takes for the field lines in the flux tube to slip pass 
through the planet by a distance equal to the planetary diameter. We consider a planet with an 
electrical conductivity a p moving relative to a magnetic field at speed v p / s = (u) p — u)*)a. If the 
planet does not drag the field at all (for example if a p = 0), the field lines would move relative to 
the planet with a linear speed v p / s . Thus, the minimum value for t max is t max = 2Rp/v p / s . On the 
other hand, if the field lines are perfectly anchored in the planet (for example when a p = oo) then 
tmax = oo. The induced field lines would wrap around the host stars with the planet's synodic 
orbit (i.e. its motion relative to the stellar spin). 

Based on extrapolation from analogous considerations (Aly 1985, Aly & Kuijpers, van Balle- 
gooijen 1994), we hypothesize that magnetic reconnection may occur when the azimuthal compo- 
nent of the induced (and "dragged" ) magnetic field well outside the planet becomes comparable 
to the unperturbed stellar dipole field. We assume the time growth timescale for the azimuthal 
component of the field to the period t syn = 2ir/(u p — w*) of the planet's synodic period. If it occurs, 
reconnection would lead to short-circuit and a burst of intense ohmic dissipation in the planet. We 
discuss the possibility of magnetic reconnection again in section (16. 4h . In the limit of infinite t max , 
the relative motion between the planet and the stellar field would restore the electric field and 
re-establish the circuit on the timescale of t& with a reduced effective conductivity (or equivalently 
an enhanced resistivity and magnetic diffusivity). We shall consider elsewhere the possibility of 
such an episodic electrodynamic process. 

In general, the conductivity a p would fall between and oo and the field is dragged without be- 
ing completely anchored on the planet. The induced field lines are distorted and partially wrapped 
around the host star. In the limit that tA < t max < t syn , it is possible to complete a steady circuit 
of unipolar induction without reconnection. 

The same argument also holds for the star with electrical conductivity dragging its own 
field lines so that the intensity of the slippage between the field lines and the planet depends on 
both the conductivity of the planet and star. Modeling the interaction between Io and Jupiter, 
Goldreich k, Lynden-Bell (1969) took into account ai a and oj U piter- Dermott (1970) included the 
contribution of the flux tube's conductivity in the expression of the slippage of the field relative to 
both Io and Jupiter. 
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In the present context, the linear speed of the slippage between the flux tube and the planet is 

Vp/s 



1 + (7£* + Ktube)/ (Tttube + TZp) 

(In the above expression, we use 7Z P instead of 21Z P as in Dermott). 
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In a complete circuit, the maximum time available for the Alfven waves to propagate from the 
planet to the star and return to the planet is 

2R P _ (1 + w) 

tmax — — 2Iip- - (32) 

Vslip (ujp-ujxja 
where w = + TZtube) / {TZp + TZtube)- It is common to neglect TZtube- 



5. Resistance and Alfven speed along the circuit 

We derived the analytical expressions of the total intensity (Section 3.4) and the time available 
for the Alfven waves to travel back and forth between the planet and the foot of the flux tube (Sec- 
tion 4.2). In order to determine their numerical values, we calculate in this section 1) the planet's 
integrated resistance TZ P , 2) the resistance TZtube and Alfven speed VA,tube along the flux tube, 3) 
and the resistance across the foot of the flux tube just below the star's surface TZ* (perpendicular 
to the magnetic field and parallel to the electric potential gradient) and the Alfven speed along the 
foot of the flux tube VA,tube- Figures [5] and [6] provide a summary of the condition of validity. 



5.1. Resistance in the Planet 7Z P 

The electrical conductivity profile of a super-Earth is unclear. We thus first discuss the (better 
characterized) conductivity profile of present day Earth. Lorrain et al. (2006) estimated the 
electric conductivity of the present-day Earth mantle to range between 10~ 2 ohm -1 m" 1 and 
10 3 ohm -1 m -1 and 10 5 ohm -1 m -1 for the inner core (also see Stevenson, 2003). Merrill et al. 
(1996, pp. 273-277) similarly argues for an electrical conductivity of about 5-8 xlO 5 ohm -1 m -1 
for the core of the Earth, and between 3 and 100 ohm -1 m _1 for the lower mantle. For the upper 
mantle (and crust), Obiekezie & Okeke (2010) calculate an electrical conductivity increasing from 
the surface (about 3 x 10~ 2 ohm -1 m -1 ) to 10 _1 ohm -1 m _1 at around 500km. The electrical 
conductivity of the Earth is therefore minimal and between 3 x 10~ 2 and 10 _1 ohm _1 m _1 ) for a 
few hundred kilometers and then increase with depth. 

In the present application, we are primarily interested in the interaction between super-Earths 
and their host stars when they are relatively young (up to a few 10 7 yr). During this stage, the 
stellar magnetic field is intense and the close-in super-Earths may be intensely heated by giant 
impacts, tidal and ohmic dissipation. Super-Earths with a molten crust are likely to have higher 
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conductivities than the present-date terrestrial planets (for example, Rikitake (1966) expressed 
the conductivity of rocks and metals on the Earth as a sum of exp(—Ei/kT), and Umemoto et 
al. (2006) estimated a conductivity at the core of a super-Earth and hot-Jupiter to be around 
10 6 ohm -1 m _1 ). The stellar radiation alone would raise the super-Earth surface temperature 
to about 1500K. Ohmic dissipation inside the planet may provide an additional source of thermal 
energy (see Section f6.3f) . In a thermal equilibrium, the planet's surface temperature may sometimes 
exceed 2000K, at which silicate melts and raise the electric conductivity to around 10 ohm -1 m _1 
(for 1400K, Waff & Weill, 1975). 

It is therefore reasonable to assume that the electrical conductivity of a super-Earth is most 
likely higher but at least that of the Earth. The electrical conductivity in a super earth would thus 
be several order of magnitude higher than 1 ohm -1 m _1 in the core and lower mantle and arguably 
also in the upper mantle. Besides, a conductivity of 0.1-1 ohm -1 m _1 (i.e. 10 times lower than the 
value we use) in an area spanning 10 % of the planet (roughly the thickness of the upper mantle) 
would at most double total the resistance. 

In addition, for a super-Earth, the characteristic speed (for example the planet linear speed 
in a frame co-rotating with the star) is much faster that for the field to diffuse across it, i.e. 
v ^ i = wnere ^ is the characteristic length. Therefore, we can neglect the diffusion (second 
term on the right hand side of the MHD equation (Equation ([T])) compared to the induction (first 
term) . 

The integrated resistance in the geometry of the planet is 

1 L 1 

Kp=q- = inr (33) 

where S is the cross section. Depending of the geometry of the current inside the planet, the 
formula could have multiplicative factors, usually of order unity. 

With these approximations, we find 

n p = —— «7x 10~ 8 Ohm (34) 
RpO 

The value we use in our fiducial calculation is !Z p = 7 x 10 _8 Ohm. If our assumption that the 
conductivity in a super-Earth is higher than that in the present-day Earth is inappropriate, this 
resistance could be higher. But it may also be much lower if the super-Earth has a substantial 
atmosphere which is extensively photo-ionized or a fully molten core where the alkali metals are 
partially ionized. 



5.2. Resistance and Alfven Speed Along the Flux Tube 



Electrodynamics along the flux tube determines the propagation of the induced electric field 
between the planet and its host star. The total resistance along the flux tube determines changes 
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in the electric potential at the foot of the flux tube. The Alfven speed determines the propagation 
speed of the disturbance. 



5.2.1. Total Resistance Along the Flux Tube 1Z tu be 

The resistance of the flux tube is also difficult to estimate accurately. Goldreich & Lynden-Bell 
(1969) simply assumed the electric conductivity to be infinite along the magnetic field lines and did 
not include TZtube i n their equations. Dermott included TZtube in the equations but, during numerical 
applications, assumed it to be negligible in front of the resistance of the satellite Io. In all previous 
investigations, conductivity across the field lines in the tenuous region between the planet and the 
star is assumed to be negligible. 

We provide here an estimate of the order of magnitude of the resistance of the flux tube. If 
we assume the plasma between the star and the planet to be fully ionized, the electric conductivity 
along (parallel) to the magnetic field line would 

c-o = — (35) 
m e v e 

where n e is the volumic number of free electrons, e and m e are the charge and mass of the electron, 
and v e the collisional frequency of the electrons with electrons and ions/protons (we assume these 
two collisional frequency to be the same). 



We take v e = — — 4^ ^ with the electron thermal speed < v e >= ■>/ 2fcgT . We thus get 

v e = n P 1.2 x our model, we consider a star with a surface 

e e u 1.2x10 °m e ' 

temperature T* = 4000K, and a planet with the equilibrium temperature T p ~ 1,500K for the 
planet with an a = 0.04 AU. For an average temperature of 2000K between the star and the planet, 
we find o"o ^ 2000 Ohm -1 m _1 so that 

L a 4.6x10-° , , 

TZtube = — g = 5 w p , ^— = 7 Ohm (36) 

(TOO {TTUp){K p Jtube)O'0 Jtube 

where jtube is between and 1 such that R p ft u be is equal to the thickness of the volumic current that 
flows along the field lines. This resistance is usually negligible compared to the other resistances 
involved in the model (especially that of the star) , except if the volumic currents are confined in an 
extremely thin layer at the surface of the flux tube. Therefore, we neglect the potential difference, 
along each field line in the flux tube between the surfaces of the planet and star. 



5.2.2. Travel Time Along the Flux Tube Between the Planet and (the Top of) the Stellar Surface 

We assume that the plasma between the star and the planet is fully ionized and estimate VA,tube 
under various different situations. 
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1) We first consider the epoch shortly after the super-Earth has migrated to the stellar proxim- 
ity through planet-disk tidal interaction. In opaque inner regions of their natal disks, super-Earths' 
type I migration generally stalls at a radius r where the has a po sitive radial gradie nt with a 



scale height Ar = Y>d/ (dYi d / dr) which is a fraction (~ 0.1 — 0.2) of r (jMasset et al 20061 ). Special 



locations include narrow transition r egions between act ive inner region and dead zone as well as 



outer edge of magnetospheric cavity (IKretke et al. 20091 ). 



We consider a super-Earth is embedded in a disk with an effective thickness Hd ~ c s /ui p ~ 
0.01 — O.lr and a steady state mass transfer rate M d = 2irY<dU a r throughout the disk where T, d 
and c s are the surface density and sound speed of the gas. Using an ad hoc a prescription for the 
effective viscosity u = aujpH'j, the radial velocity of the disk gas is U d ~ -3v/2r = -3aHjuj p /2r 
and the characteristic density at the disk mid plane is 

"^^'-edasj (3T> 

where a is the turbulent transport efficiency factor and may have an effective magnitude ~ 10~ 2 — 
10~ 3 (Hartmann et al. 1998). In untruncated protostellar disks around classical T Tauri with 
Md ~ 10~ 7 — 1O _8 M0 yr _1 , pd ~ 10~ 6 — 10~ 7 kg cm~ 3 at the edge of the magnetospheric cavity 
r ~ 0.04 AU. The corresponding Alfven speed is 

V A ,tube = 2 p ~ lOW 1 . (38) 

Near the disk inner edge, the characteristic wave propagation timescale t tu be ~ ^ r /^A,tube ~ 10 4 s 
may be too long to maintain a circuit. Note that if the super-Earth is stalled near the transition 
region between active and dead zones, tt u b e would be longer not only because this region is further 
away from the host star, but also because interior to it does not vanish. 

However, around stars with ages larger than 10 7 yr, Md may decline below that found around 
T Tauri stars and t tu be can be redu ced substantially. If the o bserved weak (or absences of) NIR 



excess around young stellar objects (jSicilia-Aguilar et al. 20061 ) in clusters with age of ~ 10 Myr is 



due to the depletion of inner holes in both gas and dust, and hence pd would be substantially 
smaller than the values estimated above. Thus in the post T Tauri phase, the Alfven speed around 
the host stars of super-Earths wound increase to sufficiently large values to enable the circuit to be 
closed, especially if we take into account the resistances in the calculation of the speed of slippage 
through the planet. 

2) Density around the flux tube does not decline indefinitely. Even after the disk is completely 
depleted or truncated in the proximity of the planet's orbit, the planet may be surrounded by a 
spherically symmetric component of stellar outflow with a speed Vfi ow and a mass flux. In this 
case, the volumic mass distribution is 

fcW -sfe (39> 
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The Alfven speed between the planet and the star at radius r (with the origin at the center of the 
star) is thus 

(40) 



V A ,tube{r) = -2 J — • 



4vr M 

Using Vfi ow = 100 km s _1 , the numerical applications give VA,tube (R* ) — 10 8 ms~ l and UA,iube(a) — 
6 x 10 6 mr 1 . The time it takes the Alfven waves to travel down the flux tube is 



Hubs 



dr — 

m 



i 



4tt M 

VOVflow 



3m \ 



4tt M 

VOVflow 



I [ — 

a 



(41) 



where the integral being for r varying from the surface of the star to the planet. If one takes 
M= lO" lo A/ yr" 1 and Vfiow — 100 km s , the numerical application then gives tA,tube — 300s. 
The magnitude of t& tube would be smaller for winds with faster speeds or lower mass loss rate. 



5.3. Resistance in the Star Across the Foot of the Flux Tube and Alfven Speed 
Along the Magnetic Field in the Star at the Foot of the Flux Tube 

The resistance of the foot of the flux tube determines the total intensity in the circuit, and 
most of the travel time of the Alfven waves occur at the foot of the flux tube. 



5.3.1. Temperature and Pressure of the Stellar Outer Layer in an Isothermal Approximation 

For the outer layer of the star, we adopt an isothermal approximation and assume a spherical 
symmetry. The pressure and temperature P(r) and T are then given by, 

T(r) = T* 
P(r) = P(R*) exp 

P(R*) = \~ 
3 K 

where k and fx are the opacity and molecular weight at the photosphere and R g = Na^b / M-H — 
8.3 x 10 3 in SI units, with Ma the Avogadro number, ks the Boltzmann constant, and MIh the 
molar mass of the hydrogen atom. The volumic mass can also be calculated using the ideal gas 
equation of state. For the models presented here, we neglect any change in T(r) and P(r) due to 
the local ohmic heating at the foot of the flux tube. Discussions in Section 16.21 shows the possible 
existence of a hot spot at the foot of the flux tube. Self consistent treatment of a potential feedback 
effect will be analyzed elsewhere. 



GM*n (R* 
R g T(r)R, V r 



(42) 
(43) 

(44) 
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5.3.2. Conductivity in the Stellar Atmosphere 

The details of the derivation of the conductivity in the stellar atmosphere (foot of the flux 
tube) are given in Appendix A. In Figured! we show that, at the foot of the flux tube, current flow 
across the field lines, as a series of parallel circuits. We calculate the effective resistance and Alfven 
travel timescale to determine the depth of penetration. We use the Saha's equation to derive the 
ionization fraction. Following Fejer (1965) we refer ao as the electric conductivity parallel to the 
magnetic field lines, and a p = 1+ ^J°/ V p the Pedersen conductivity parallel to the electric field. We 
define oo e = eB/m e to be the electron gyro- frequency, v e to be the mean collision frequency of the 
electrons with the neutral gas (see Equations (A7) and (A8)) and r= to be the radius at which 

u e (r=) = w e (r=). (45) 

We find r= ~ 1.3962 x 10 9 m (given with several significant figures as an intermediate value in the 
series of numerical applications). Since co e (r) « v e (r~) at r < r= and u e (r) > v e { r ) at r > r = , we 
write the Pedersen conductivity 

a p {r < r = ) = a (r) (46) 
a p (r>r = ) = <r (r)(^) ■ (47) 

In contrast to the region between the planet and its host star, gas in the stellar atmosphere is 
partially ionized. Substituting the appropriate value for o"o from Equations (A4) and (A5) we find 

a ^ 2 = m exp (™m) r ° iP(r)f ' 2 m 

hi - in- 1 ' (t2.Y 2 (2irm«) 3/4 -i/4 (IJ&m \ 1/2 
Ql 



w ~"{fJ ^W^ k ° {^ (50) 



where the numerical values of the constants in SI units are L a = 6.17 x 10 6 , Q a = 2.93 x 10" 4 , and 
L a /Q 2 a = 7.2 x 10 13 . 



5.3.3. Alfven Speed and Resistance 
In Equation (jlip . we showed that is given by 

^* = (51) 
2hs 

where s has been previously defined as cosOp. Here, we decompose S, the integral from r pn to 
i?* of the electric conductivity, into two parts Si and £2, respectively the integral of the electric 
conductivity a p (z) from r pn to r= and from r= to i?* such that 
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(jp(r)dr 



a p (r)dr = L 



T ( r )3/4 



-.exp 



E 



L a 



<jp(r)dr = -p^exp 



Ql 



E 



2k B T(r 



where y\ = {-¥■) (R*/a)' i/z and 1/2 = 2R p {R*/a) ilz are respectively the length and width of the foot 
of the flux tube, r= is defined above, and the penetration radius r pn is to be determined below (see 
Figure^]). The numerical application gives S2 ~ 1-5 x 10 4 ohm -1 , and the analytical expression of 
Si depends on r pn . Nevertheless, in our fiducial model, the numerical value of Si does not depend 
significantly on r pn and we get Si w 4.3 x 10 4 ohm _1 . These values lead to = 8.6 x 1(T 6 ohm. 

At the foot of the flux tube on the surface of the star, the volumic mass is g(r) = P(r)m p /k B T 
and the expression for the ionization x is given in the appendix A. Therefore, the Alfven speed in 
the stellar atmosphere and the time t& is thus 
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2k B T(r)j m 2 T{r) l l 4 
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dr (53) 
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(55) 

(56) 

(57) 
(58) 



In SI units, W = 0.038, B(r) = 2cos9 Florin/ {An R^) (with cosQf ~ 1, and the integral going 
from r pn to R* ( Vp n is the radius that determines the effective height R* — r pn of the foot of the flux 
tube). Using the fiducial values for the parameters, we get va,» — 2.6 x 10 8 exp [-717.5(^-1)]. 
Clearly, va,* decreases sharply from the surface toward the interior. r pn is thus the smallest radius 
that still enables the model to be valid (i.e. the deepest that a perturbation of the field line can 
penetrate inside the star and back to the planet in less than t max (see Figure [8]) . 



The time it takes for the Alfven wave to travel from the surface of the star to the bottom of 
the flux tube is 

f dz , . 

U* = / — (59) 



Jz VA,*(z) 

We then equate the total travel time 2(tA,* + ttube) (there is a coefficient "2" since the wave 
needs to go from the planet to the star and back to the planet) defined in Equations (|59p and (|4ip 
with the total time available t max defined in Equation (|32p 



2(tA,* + ttube) — t f 



(60) 
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and solve for r pn . The numerical application gives r pn ps 1.3718 x 10 9 m for a fast rotating star and 
r pn ~ 1.3746 x 10 9 m for a slow rotating star (given here with several significant digits simply as 
an intermediate value in the thread of numerical applications). Having determined r pn , one could 
now calculate self-consistently Si and 7£*. 

6. Numerical Applications and Discussion 

The quantities we have determined above are applicable for the fiducial model we adopted 
here. A more general application (for host stars of different masses) will be presented elsewhere. 

6.1. Numerical Applications 

The numerical values of the intensity, ohmic dissipation, and torque in the star and the planet 
respectively for a slow rotating star (with a spin period 8 days) and a fast rotating star (with a 
spin period 0.8 days) are given below. The values below for U and I correspond to one of the two 
circuits (each circuit has the same value of U and I) , and the values for V and T are for the entire 
planet and entire star (both circuits combined). 

U = 6.7 x 10 9 V and 2.8 x 10 10 V (for slow and fast rotator respectively) 

/ = 7.8 x 10 14 A and 3.2 x 10 15 A 

7?„ = 10 25 W and 2 x 10 26 W 

V p = 4x 10 22 W and 7 x 10 23 W 

| % |=| T p |= 6 x 10 29 N m and 2 x 10 30 N m. 

6.2. Ohmic Dissipation at the Foot of the Flux Tube on the Star and Hot-Spots 

A super-Earth orbiting at a = 0.04 AU from its host star induces an ohmic dissipation at the 
foot of the flux tube on the surface of the star of 5 x 10 24 W (for a stellar spin period of 8 days) 
and 9 x 10 25 W (for a spin period of 0.8 days, which is about 1 to 5% of the stellar luminosity 
(L* : total = 3.6 x 10 26 W). This effect would create an observable hot-spot at the surface of the 
star. The rate of energy dissipation per surface area at the foot of the flux tube would be about 
2 x 10 10 W m~ 2 which is three orders of magnitude higher than the intrinsic radiative flux from the 
surface of a typical T Tauri star. 

Since Si > S2, most of the dissipation occurs in the region between r pn and r = . We note 
that the density scale at the photosphere (-R*), 5r* = RgT^R^/GM^fi ~ 5 x 10 5 m is much smaller 
than - r= = 3.8 x 10 6 m and R* 

~ r pn — 1.55 x 10^m. When the dissipated energy emerges 
from the stellar photosphere, the actual area of the hot spot may be diffused to several times the 
area of the foot of the flux tube. The corresponding temperature of the hot spot would be ~ 2 — 3 
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that elsewhere on the stellar surface. In our model, we consider a star with a surface temperature 
T* =4000K. Ohmic dissipation at the foot the flux tube increases the local ionization, conductivity, 
current, and torque. We will construct a self consistent model in a follow-up analysis. 



6.3. Ohmic Dissipation in the Planet and the Induced Mass Loss 

In Paper I, we considered the structural adjustment due to ohmic dissipation. In this paper, we 
have not yet considered the structure adjustment of super-Earth structure due to ohmic dissipation. 
For the slow-rotator model, the rate of ohmic dissipation is six times that the planet receives from 
its host star's irradiation. In the absence of any structural adjustment, the super-Earth may attain 
a thermal equilibrium with an effective blackbody temperature T p ~2,300 K (and much higher 
for the fast-rotator model). With this temperature, planet's core crust would surrounded by an 
ocean and an extensive atmosphere where water and hydrogen molecules readily dissociate but 
their ionization fraction remains negligible. The local density scale height is 5r p = XR p where 

A = RgTpR p /GMpfj, ~ 1/17^. (61) 



The density scale height of the hydrogen atoms is much larger than that of all other elements 
including carbon, oxygen, and silicates. The mean free path for a hydrogen atom to collide with an 
heavy element with a density n z is Ih—z = l/{n z A) where A ~ 1CP 19 m 2 is a typical cross section. 
Within 5rde ~ 2 — 3 hydrogen atoms' scale heights there are so few heavy-elemental atoms left, 
they essentially become thermally decoupled, i.e. Ih—z > Sr^e for hydrogen atoms. We refer this 
location as the decoupling radius r^ e = R p + Sr^e ~ (1.1 — 1.2) R p and the local hydrogen density 
at r^ e as Pde)- The local temperature T de = T(r^ e ) is set by the blackbody temperature T p of the 
heavy elements in the planet's photosphere. The magnitude of pd e = p(fde) for a rocky or icy planet 
can be estimated to be 

Pde = 3/tfM p /47Ti^exp - 5 de ~ 10 kg m" 3 (62) 

where typical fractional abundance of the hydrogen atoms fn < 0.1. 

Next, we consider the possibility of significant loss of hydrogen atmosphere. Planetary outflow 
is usually analyzed in the limit that atmosphere is heated by stellar irradiation. For the present 
configuration, simple estimates indicate that the hydrogen atmosphere is opaque to incident ion- 
izing photons from the host star, i.e. they are mostly absorbed by hydrogen atoms in the upper 
atmosphere. Provided hydrogen atmosphere remains mostly atomic, most visual stellar photons 
would stream pass it and be absorbed by heavy elements near R p . Transit light curves of such a 
super Earth i n Ly a photons would be m uch deeper than that for visual photons, as in the case of 



HD 209458b, (IVidal-Madjar et al. 20031 ). Thus, most regions of the atmosphere is not affected by 



either ohmic dissipation or irradiation. 



The most important heat input is the ohmic dissipation which takes place at the base of the 
atmosphere. If we neglect energy deposition and loss in the atmosphere, our problem would reduce 
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to a simple spherical Bondi (or Parker) solution. In a steady state, the governing continuity and 
momentum equations at a location r would reduce to 



M H = Airp H V H r 2 (63) 
T , dV H c 2 s d PH GM p 3GM*r 

V H-^- = o — + 3 ( 64 ) 

or pn or r z or 



where the last term represents host star's tidal force (jDobbs-Dixon et al. 20071 ) and Vh is the radial 
velocity. Together they reduce to 

(Vg - c 2 s ) dhiV H _ 2c 2 s GM p 3GM*r 

en ~~ 5 ' ^ • V"^J 

r omr r r z a 6 



Transonic point (where Vh = c s ) occurs (|Lubow fc Shu 197a ; iGu et al. 20031 ) near the Roche radius 
r ~ Rr = (Mp/3M*) 1/ ' 3 a. Interior to the transonic point, flow is subsonic. 

In order to make further progress, we need to estimate the energy budget of the atmosphere. At 
the base of the atmosphere, ohmic dissipation occurs primarily due to collision of charged (provided 
by the heavy elements) and neutral particles. Most of the dissipated energy is emitted to space at 
R p as blackbody radiation. Below r^ e , hydrogen atoms attain T p (~ 2,300K) through conduction 
as all other particles. Above r^ e , hydrogen atoms attain different density distribution. 

For computation simplicity, let us first consider an isothermal equation of state. For an analytic 
approximation, we neglect the advection contribution to the momentum equation and obtain 

p(r) ~ p de exp(r de /Xr - 1/A). (66) 

At large distances (r > (2 — 3)r^ e ~ 3i? p ) but still well within Rr(~ 10 10 cm ~ 7R P ), hydrogen's 
density app roaches to an i sochoric limiting value p^ ~4x 10~ 8 pd e ~ 10 -7 kg m~ 3 . The mass loss 



rate at Rr (|Li et al. 2010] ) becomes 

Mhydro - ^RrPooCs ~ 5 x 10 12 kgs -1 . (67) 
At this rate, the total hydrogen mass in the planet JhM p is depleted in a few Myr. 

Note that hydrogen atoms contained within this region is negligible compared with M p as 
we have assumed in the momentum equation. In addition, the collisional mean free path be- 
tween hydrogen atoms Ih-h = mH/(pA) is small compared the density scale height 5r p and 
more importantly Rr. In this limits, it is more appropriate consider outflow in the hydrodynamic 



limit (IMurrav-Clav et al. 200911. a s we have done above, rather than use the Jeans' escape formula 



(jLecavelier des Etangs et al. 20061 ). 



If the planet's atmosphere is maintained above the recombination temperature so that it is 
primarily compose of hydrogen atoms, the main cooling process would be the emission of Ly a 
photons at a rate 

A ~ 7.5 x 10 _20 xn? r exp - (1.2 x 10 5 K/T)Jm~ V 1 . 
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After integrating over the entire volume ~ 4-7ri?|j/3, the total energy loss rate is LL ya ~ 10 16 x Watt 
which is substantially below fraction of the dissipated energy flux (fn'Pp) carried by the hydrogen 
atoms. (In the above estimate, we use the asymptotic value of p^ to estimate A.) However, with 



advected with the escape hydrogen gas. 

Based on hydrogen atoms' ineffective absorption and emission rates, it is natural to contem- 
plate the possibility that the planet's atmosphere expands adiabatically. In the limit that ohmic 
dissipation provides the only source of heating at its base, planet's atmosphere may be convectively 
unstable. Efficient convection also leads to constant entropy. 

Using the conventional polytropic approximation (in which P = Kp' 1 and 7 = 1.4) for an 
adiabatic hydrogen atmosphere, a stationary quasi hydrostatic solution can be constructed with 



The above equation implies that with an adiabatic equation of state, both density and temperature 
in the hydrogen atmosphere vanishes within Sr^e/ij — 1) ~ 2.5<5rrf e . Unless planet's photospheric 
radius can expand (see paper I) significantly, there would be no outflow, despite the intense ohmic 
dissipation below R p , and all the thermal energy generated would efficiently radiated by atomic 
emission from heavy elements. 

However, planet's atmosphere may be heated to prevent its temperature from plummeting 
below that (~ 2000K) for hydrogen molecules to recombine. Rotational and vibrational bands of 
hydrogen molecules not only provide emission mechanisms but also opacity sources to absorb the 
incident stellar irradiation and to diffuse thermal energy in the planet's atmosphere from its heated 
base to its upper layers. In the mildly heated case, we anticipate the planet's hydrogen atmosphere 
to attain an equilibrium temperature so that the incident deposition of stellar photon energy would 
be balanced by planet's reprocessed luminosity. (For our fiducial model, the equilibrium temper- 
ature is T p ~ 2000K.) In this limit, the loss of planet's hydrogen atmosphere relies more critically 
on the atmosphere ability to maintain a shallow temperature gradient than that to generate energy 
through ohmic dissipation . This situation has already been analyzed in the context of HD 209458b 
dMurrav-Clav et al. 2009T ). 

However, enhanced sources of energy may also expand the radius of planet's photosphere 
well beyond R p . This is a distinctive possibility for the fast-rotator model in which case P p is 
another 18 times larger. This increase in the ohmic dissipation rate is due to the relatively large 
differential motion between the planet and the magnetosphere of its host star. If the planet's 
photosphere remains at R p , the enhanced energy source would increase T p by a factor of 2 which 
is comparable to the magnitude of T*. At this temperature, opacity due to H~ process becomes 
significant. Planet's envelope and photosphere may well expand, leading to a possible runaway 
ohmic heating. The above discussion clearly warrants further discussions and detailed treatments 
of radiation transfer in this type of super Earths We shall carry out and present these analysis in 



the magnitude of Mhydro in Equation ([6T|h we find that a significant fraction of fn'Pp may be 




(69) 
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a future paper. 

The loss of planet's atmospheric hydrogen is likely to occur on a more rapid pace. It remains 
to be demonstrated that for the intense heating cases, how far up in the atmosphere does thermal 
decoupling between hydrogen and heavy elements occur. If the planet's photosphere is well within 
Rr, the density scale height of most other heavy elements such as carbon and oxygen above are 
sufficiently small that they may be effectively retained near R p . Oxygen atoms may combine with 
Mg, Fe, Ca, Na, Al, and Ti silicates to form high-density minerals such as enstatites, olivines, and 
pyroxenes. Planets composed mostly such substances are expected to have compact sizes (Valencia 
et al. 2010). Thus, it is likely that super Earths which migrated early to the proximity of their 
strongly magnetiz ed host stars may attain relatively compact sizes as in the case of COROT 7-b 



(jLeger et al. 20091 ) and planets around Kepler 11. 



The rate of ohmic dissipation in short-period super Earths is likely to diminish as their host 
stars magnetic field weakens with age. As their semimajor axis increases, planets which undergo 
outward migration around rapidly spinning host stars also encounter less intense stellar dipole 
field. Some residual oxygen atoms in the atmosphere may recombine to form oxygen molecules 
during the decline of the ohmic dissipation rate. Oxygen molecules are particularly important 
because they have been suggested as a bio-marker for the detection of life elsewhere in the Universe 



(jDesMarais et al. 2002 ). 



6.4. Discussion about Some Approximations 

We presented here a preliminary model for the unipolar induction model. Some approxima- 
tions were made for computational convenience, and we briefly discuss here the validity of the 
approximations which have not yet been discussed in the paper. 

We only took into account hydrogen for the calculation of the conductivity/resistance of the 
star. In a realistic model, especially for low-mass stars, other elements may become important 
contributors of the ionization fraction. A more comprehensive study will be presented elsewhere. In 
addition, we suggested that the foot of the flux tube at the stellar atmosphere would be significantly 
heated. However, we used T=4000K for the stellar surface temperature (usual T Tauri star). The 
feedback on the stellar temperature due to the circuit may be included in later models. 

For the evaluation of the planet's electrical resistivity. The value we adopted (7 x 10~ 8 ohm) 
seems to be the most uncertain value in our calculation. This value of the resistance of the planet 
used here is likely to be a lower boundary for the mantle of the planet (a metallic core might have 
even higher conductivity). If the real resistance where to be lower, then 1) the time t max available 
for the Alfven waves to travel around the circuit would increase, which would result in a deeper 
foot of the flux tube and would also enable the model to hold for larger semi-major axes and 2) 
the ohmic dissipation in the planet would decrease. Nevertheless, an increase in the depth of the 
foot of the flux tube would not affect much the total resistance of the foot of the flux tube 1Z* . 
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Induction at the foot of the flux tube. Since the conductivity is very high along the magnetic 
field lines, the plasma in the star's magnetosphere rotates with the magnetic field lines. Therefore, 
the plasma contained in the flux tube also moves with the magnetic field lines as they are dragged 
along by the planet, and thus moves relative to the unperturbed magnetic field line. Therefore, 
just as the induction in the planet is due to the relative motion between the frame co-moving 
with the planet and the frame rotating with the magnetosphere, there can also be a magnetic 
induction in the plasma enclosed by the foot of the flux tube. The order of magnitude of this 
phenomenon will be at most comparable to the order of magnitude of the phenomenon presently 
described. Ferraro & Plumpton (1966) provide a brief discussion of the problems raised by two 
good concentric conductors rotating in a magnetic field at different angular speeds. 

As mentioned earlier, there may be magnetic reconnection if the induced field dominate over 
the unperturbed stellar dipole field. The field lines also tend to wrap around the planet when 
the synodic period (T syno( i = - ) is small compared to t max , the time required for the field 
lines constituting the flux tube to move across the diameter of the planet. Using (|32|) . we find 
that t max = + w)T syno d with w = ^*+tiTT\ ~ TIT' ^ or our P arame ters, this corresponds 

to t max ~ 0.091T s?/no rf (using 1Z P = 7 x 1CT 8 ohm according to (f3~i|) and 1Z* = 8.6 x 1CT 6 ohm 
according to section (|5.3.3j) ) . 

This also provides an upper limit on the ratio of resistances w in order to stay with a model 
without reconnection. Indeed, magnetic reconnection may occur when t max is larger than a few 
T synod, an d one would thus arguably stay in the regime without frequent magnetic reconnection 
when 

ira 

l + w< K— (70) 
Rp 

with K larger than 1, and w as defined above. Using R p = 2R§ and a= 0.04AU, we find ~ 1300. 
We have neglected the resistance of the flux tube Tttube i n front of 1Z P and 1Z* but this approximation 
may break down in extreme cases. Nevertheless, in most cases, Equation ()70|) means that 1Z*/1Z P is 
smaller than 1300. The resistance of the planet depends mainly on its composition, and structure, 
which would adjust to the strong ohmic dissipation in its interior. The resistance of the foot of the 
flux tube in the stellar atmosphere would depend on the metallicity and the temperature, which 
would also adjust to the strong ohmic dissipation. 

Previously in the paper, we also discussed that the travel time of the Alfven wave can be at 
most t max (Equation ([60)) ), 7£*, through its relationship with the variable depth of penetration, 
would self-consistently adjust depending on the parameters of the model. Indeed, larger 1Z* leads 
to larger w, then larger t max , and thus larger depth of penetration since the Alfven waves have 
more time to travel between the planet and the star along the flux tube, into the stellar atmosphere 
at the foot of the flux tube, and back to the planet. Deeper depth of penetration then results in 
smaller 1Z* (equivalent resistance with resistances in parallel). 

The value 7Z P is less directly constrained by the model although it of course depends on 
the parameters chosen for the model. Nevertheless, changes in the value of 1Z P would result in 
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adjustments in 7Z* through the mechanism mentioned just above. 

Goldreich and Lynden-Bell also interpret the torque calculated above (section 16. ip in terms of 
a toroidal magnetic stress due to a distortion of B^ in the azimuthal direction, i.e. the direction of 
the motion. Neglecting the induced field in the r — z (meridional) direction, they determined the 
longitude of the flux tube from the ratio of the induced and the unperturbed stellar-dipole field. 
They then determined, for the Jupiter-Io system, the forward-sweeping angle (or the backward- 
sweeping angle in the case of a slowly rotating star) of each field line as it leaves the Io to be 13 
deg. When a similar approach is adopted in the present model, we find this angle may be close 
to 90 deg. This large distortion is due to a strong torque induced by the unipolar circuit with a 
relatively small 1Z* (and thus a large intensity). For such a large field distortion, Goldreich and 
Lynden-Bell suggested that the induction circuit may be broken by field reconnection. We shall 
further examine this possibility elsewhere and determine whether it may significantly weaken the 
effective torque. 



7. Summary and discussions 

With the advent of high-precision radial velocity and transit surveys, we have entered an era of 
super-Earth discovery. Although the detection probability (due to observational selection effects) 
decreases with planets' period, three times more planets are found with period between 3 and 10 
days than between 1 and 3 days. We suggest that super-Earths' interaction with the magnetosphere 
of their host stars may be one possible mechanism for this dichotomy. 

In this paper, we analyze the electrodynamics of super Earths orbiting in the proximity of 
strongly magnetized T Tauri stars. We constructed a fiducial model in which the planet's orbital 
frequency is not synchronized with the star's spin. Their relative motion enables the planet to 
continually encounter field lines which are locked on the star. As a good (but not perfect) conductor, 
an emf is induced across the planet (along the semi-major axis). We estimate planet's conductivity 
and show that the stellar fields slip through the planet with a drift speed considerably slower than 
its Keplerian speed. 

We show that conductivity along the field line is likely to be large and the perturbed potential 
(due to the induced electric field) propagates along a flux tube away from the planet with an Alfven 
speed. We show that for planets with period less than 3 or so days, the disturbance can reach the 
surface of the star and return before stellar fields have drifted through the planet. 

The foot of the flux tube is implanted to stellar surface. As density increases with depth below 
the photosphere, the Alfven speed decreases. Penetration depth of the flux tube is determined by 
the condition that the timescale required for Alfven waves to complete a circuit between the planet 
and its host star is comparable to that for the stellar field to drift through the planet. 

Across the foot of the flux tube on the stellar surface, the potential drop induces a current 
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to flow across it. We show that the resistance on the surface of the star is larger than that in the 
planet. Consequently, the intensity of the current is determined by the resistivity on the star. We 
quantitatively determine this resistivity, the associated current, ohmic dissipation rate, and torque 
due to the Lorentz force. The ohmic dissipation in the star at the foot of the flux tube could also 
induce an observable hot-spot. 

The source of energy is the differential motion between the planet and the magnetosphere of 
its host star. The Lorentz force on the planet and its host star leads to an evolution toward a state 
of synchronous rotation. Inside the corotation radius, planets tend to lose angular momentum and 
migrate inward and the opposite trend occurs outside the corotation radius. Consequently planets 
inside corotation migrate inward and those outside corotation migrate outward. 

For super-Earths with periods less than 3 days, the timescale for orbital evolution can be 
comparable or shorter than a few Myr (the timescale over which intense stellar magnetic field is 
maintained). The low abundance of super-Earths with period less than 3 days may be due to their 
infant mortality. 

Due to their finite conductivity, ohmic dissipation also occurs on within the super-Earths. The 
heating rate depends on planet's poorly determined resistivity. Its magnitude can be comparable 
to or larger than that the planet received from the stellar irradiation. The intense rate of ohmic 
dissipation may cause water and hydrogen molecule to dissociated and hydrogen atoms to segregate 
from other heavy elements. It is unclear whether a substantial fraction of the hydrogen atom may 
escape though hydrodynamic outflows. As the field decay with maturing stars, remaining excess 
oxygen atoms may either be incorporated in high density minerals or form oxygen molecules. Either 
of these processes can lead to consequences which may be observable in the near future. 

There are several uncertainties which warrant further investigation. Conductivity in super- 
Earths and their host star need further study. We have not yet apply these results to a wide range 
of stellar and planetary models. The effect of feedback due to the adjustment of planet's and star's 
heated atmosphere also need to be examine. Perhaps the largest uncertainty is whether the intense 
induced field can lead to magnetic reconnection and the breaking of the circuit. Reconnection 
would increase the effective magnetic diffusivity and severely weakens the effective torque. 

In order to directly compare with observations, we also need to consider a diverse range of 
planetary orbits. For example, this process may not work for planets with period longer than a few 
days. Finally, it would be of interest to determine whether the intense electromagnetic interaction 
between super Earths and their host stars can be directly observed in the radio-wave frequency 
range. 

Nevertheless, we show that electro dynamic interaction is an important process for the orbital 
and structure evolution of super-Earths as well as hot Jupiters. Along with many other physical 
processes it introduces diversity in the present-day configuration of extra solar planetary systems. 
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A. Electric conductivity at the foot of the flux tube 

The Saha's equation gives the ionization fraction of the hydrogen atom x = \J i+k h » wnere 

1 (2>rm e )3/2 5/2 __( E\ 



K »^=W)-i^ {m exp {-w) (A1) 

. where P(r) is the pressure, m e the electron mass, h Planck's constant, and E the ionization 
energy of the hydrogen atom. Kh is a function of r which, in the isothermal region, decreases as 
one moves from the surface of the star toward the interior. Since, in the situations considered in 
this paper, Kh is small compared to unity at the surface of the star, we get the following expression 
for the ionization rate x everywhere in the isothermal region 



x(r) « y/K H (r). (A2) 

Using the formulas given by Fejer (1965), we calculate the electric conductivity profile in the stellar 
(isothermal) outer layer. The conductivity cro( r )> which determines the current parallel to the 
magnetic lines of force, is given by 



iq e 2 / 97rm e . , , 

, oW=lW 10--^^ (A3) 

Using (1A1I) and (1A2I) . we obtain the following expression for o"o 

/ —E \ T(r) 3 / 4 
a M = L a exp ' (A4) 

V ' 1 \2kT{r)) Jp(rj V ' 

L - 1Q 19 e 2 / 97rm e '(27rm e ) 3 / 4 5/4 , . 

which decreases as one moves from the surface of the star toward the interior. The numerical value 
of the constant L a in SI units is L a ~ 6.17 x 10 6 

The (Pedersen) conductivity cr p (r), which determines the current parallel to the electric field, 
is given by 

<J P = (A6) 

where uj e (the gyro-frequency of the electron) and v e (in the limit of a gas with low ionization 
fraction, z^ e , is related to the mean collisional frequencies of the electrons with molecules of the 
neutral gas, see Draine et al. 1983) are given by 

Mr) = ^ (A7) 

vAr) _ KT" „f 128M '' r ))' / ' = 10-'f(r) ( 12 l, X' 2 m 
y ' \ 97rm e J V ; \97rm e kT(r) J v ' 
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with n the number density of neutral particles. Since the ionization rate is small, n is also the 
number density of particles which is equal to P/kT for a perfect gas. Using (|A7|) and the expression 
of a dipole magnetic field B s (r) = fiQ m/(4-7rr 3 ) (with m being the stellar magnetic moment), we 
obtain 

v e r 3 P(r) 

^ e Un iq f 9irm e kB \ 

Q " = ss'(iir) < A10 » 

where the numerical value of the constant Q a in SI units is Q a = 2.93 x 10 -4 . 
In order to compare ^(r) with unity, we define r= such that 

^(r=) = 1 (AH) 

or, equivalently, 

P(r = ) = 9m ^ y . (A12) 

Using the numerical values for the star listed above, we deduce r= = 1.3962 x 10 9 (note that in 
paper I, what we defined r= to be the transition between the isothermal and the polytropic region 
in the star, which is unrelated quantity defined here). 

Since P(r) increases rapidly when r decreases (from the stellar surface inward), one may 
distinguish two regimes by 

a p {r<r = ) = a (r) (A13) 

a P (r>r = ) = ^ (A14) 

( 

i.e. 

/ —E \ T(r) 3 / 4 

^ rSr -' - L ' exp [^sm) 7m <A15> 

Mr > r.) = | e* (^) r« (P W )»/* (A16) 

where the numerical value of the constant in SI units is L a jQ\ = 7.2 x 10 13 . 



B. Another estimate of the Alfven speed along the flux tube 



The volumic current passing through the tube is J = ev e f g/m p . In this expression, we adopt 
the current propagation v e to be the thermal speed of the electrons m e v1/2 = k^T and we define / 
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to be the fraction of gas particles that are ionized. The total flux J is the total intensity I divided 
by the cross section of the flux tube through which the current passes. It is a fraction g 2 of the 
cross section ttR^ of the planet. With these notations, we get J = I /g^R^. 

We thus obtain the following Alfven speed and tA,tube 



VA,tube 



tA,tube 



Mom 

47T 



1 vri?2 e /2k B T\ 



Ho I m p \m e J 



2 / irlm p ( m e \ 
mV ^Rle \2k B T) 



1/4 








1 - 











r J92 



which gives, for fg 2 = l, t A>tube ~ 450s. 
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Scenario 3: If the Alvfen wave propagates deeper than d pn 




Fig. 8. — Penetration depth of the unipolar induction circuit. The propagation of the induced 
disturbances from the planet to the stellar surface is tA,tube- Below the stellar surface, gas density 
increases exponentially and the Alfven speed decreases accordingly. The penetration depth of the 
planet's induced field is determined by the requirement that the timescale for the Alfven waves to 
complete the circuit equals to the field slippage timescale across the planet. The net resistivity at 
the stellar surface is determined by the gas across the flux tube in the form of parallel resistors. 



